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When  a  ball  with  a  concentric  hole  through  it  is  suspended  within  a  spin¬ 
ning  tubular  projectile  the  ball  will  rotate  in  such  a  way  that  its  hole  is 
aligned  with  the  axis  of  spir  of  the  projectile.  The  ball  thus  serves  as  an 
automatic  obturator--plugging  the  projectile  while  the  two  bodies  are  within 
the  gun  barrel  and  rotating  to  open  the  tubular  passage  following  exit  from 
the  muzzle.  The  motion  of  the  ball  is  gyrodynamic  in  nature  and  highly  depen¬ 
dent  upon  the  external  moments  on  the  ball  that  arise  because  of  its  motion'-- 
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relative  to  the  spinning  projectile. 

This  report  presents  the  results  of  an  analytical  and  experimental 
study  to  provide  an  understanding  of  such  motions  in  as  complete  a  sense  as 
is  possible.  The  general  equations  of  motion  are  formulated  together  with 
appropriate  mathematical  expressions  for  the  external  moments.  These  moments 
have  their  origins  in  the  forces  acting  on  the  ball  (gravity  and  aerodynamic 
forces  are  considered)  which,  in  the  presence  of  relative  motion  between 
the  ball  and  the  projectile,  lead  to  sliding  friction  torques.  Torques 
due  to  fluid  shear  are  also  considered  but  these  are  found  to  be  relative¬ 
ly  insignificant. 

Exact  solutions  are  obtained  using  standard  numerical  techniques.  In 
addition,  a  linear  form  is  developed  and  these  solutions  lead  to  useful 
approximations  that  are  valid  over  a  broad  range  of  operating  conditions. 

Experiments  are  described  in  which  the  general  validity  of  the  theoreti¬ 
cal  models  (exact  and  linear)  is  demonstrated  for  cases  in  which  there  is 
little  or  no  aerodynamic  load.  The  experiments,  together  with  the  linear 
approximation,  lead  to  a  semi-empirical  method  for  determining  the  effective 
coefficient  of  sliding  friction  for  such  systems. 

A  model  is  proposed  to  account  for  aerodynamic  loading  and  sensitivity 
studies  are  conducted  to  determint  the  nature  and  scope  of  the  influence 
of  various  design  parameters  upon  systme  performance.  The  linear  approxi¬ 
mation  is  proposed  as  a  useful  design  guide  when  applied  with  due  awareness 
of  its  limitations.  In  addition,  a  design  criterion  is  presented  by  means 
of  which  it  is  possible  to  avoid  designs  that  lead  to  operation  in  a 
"hovering"  region.  Such  operation  leads  to  prolonged  delays  in  ball  opening 
time.  Application  of  the  design  criterion  leads  to  ball /projectile 
designs  that  are  optimum  from  the  point  of  view  of  minimum  ball  opening  time. 
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When  a  ball  with  a  concentric  hole  through  it  is  suspended  within 
a  spinning  spherical  cavity  the  ball  will  rotate  in  such  a  way  that  its 
hole  is  aligned  with  the  axis  of  spin  of  the  cavity.  If  it  is  mounted 
within  a  spinning  tubular  projectile,  the  ball  can  serve  as  an  automatic 
obturator-plugging  the  projectile  while  the  two  bodies  are  within  the 
gun  barrel  and  rotating  to  open  the  tubular  passage  following  exit  from 
the  muzzle.  The  motion  of  the  ball  is  gyrodynamic  in  nature  and  highly 
dependent  upon  the  external  moments  on  the  ball  that  arise  because  of  its 
motion  relative  to  the  spinning  cavity. 

This  report  presents  the  results  of  an  analytical  and  experimental 
study  to  provide  an  understanding  of  such  motions  in  as  complete  a  sense 
as  is  possible.  The  general  equations  of  motion  are  formulated  together 
with  appropriate  mathematical  expressions  for  the  external  moments.  These 
moments  have  their  origins  in  the  forces  acting  on  the  ball  (gravity  and 
aerodynamic  forces  are  considered)  which,  in  the  presence  of  relative 
motion  between  the  ball  and  the  projectile,  lead  to  sliding  friction  torques. 
Torques  due  to  fluid  shear  are  also  considered  but  these  are  found  to  be 
relatively  insignificant. 

Exact  solutions  are  obtained  using  standard  numerical  techniques. 

In  addition,  a  linear  form  is  developed  and  these  solutions  lead  to  use¬ 
ful  approximations  that  are  valid  over  a  broad  range  of  operating  condi- 
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Experiments  are  described  in  which  the  general  validity  of  the 
theoretical  models  (exact  and  approximate)  is  demonstrated  for  cases  in 
which  there  is  little  or  no  aerodynamic  load.  The  experiments,  together 
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with  the  linear  approximation,  lead  to  a  semi -empirical  method  for  deter¬ 
mining  the  effective  coefficient  of  sliding  friction  for  such  systems. 

A  model  is  proposed  to  account  for  aerodynamic  loading  and  sensi¬ 
tivity  studies  are  conducted  to  determine  the  nature  and  scope  of  the  in¬ 
fluence  of  various  design  parameters  upon  system  performance.  The  linear 
approximation  is  proposed  as  a  useful  design  guide  when  applied  with  due 
awareness  of  its  limitations.  In  addition,  a  design  criterion  is  presented 
by  means  of  which  it  is  possible  to  avoid  designs  that  lead  to  operation 
in  a  "hovering"  region.  Such  operation  leads  to  prolonged  delays  in 
ball  opening  time.  Application  of  the  design  criterion  leads  to  a 
ball/projectile  design  that  is  optimum  with  respect  to  ball  opening  time. 
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INTRODUCTION 


BACKGROUND 

As  early  as  1858,  the  potential  benefits  of  tubular  projectiles  were 
recognized  for  specific  purposes.  In  that  year,  Joseph  Whitworth  (Whitworth 
Threads)  wrote  about  and  included  an  illustration  of  the  tubular  projectile 
in  the  section  on  Rifled  Firearms  of  his  Miscellaneous  Papers  on  Mechanical 
Subjects  [1] 1 .  The  projectile  devised  by  Whitworth  was  hexagonal,  with  a 
circular  hole,  and  he  noted  its  particular  "...effectiveness  in  perforating 
elastic  materials  which  prevented  them  from  closing  up."  Whitworth  further 
observed  that  the  tubular  projectile,  which  utilized  a  wooden  sabot,  penetrated 
deeper  into  masonry  than  any  with  which  he  was  acquainted. 

According  to  Charters  and  Thomas  [2],  the  "Krnka-Hebler"  projectile 
was  reported  in  the  Allgemeine  Schwei tzerische  Militarzeitung  as  having  been 
so  successful  that  the  Ordnance  Department  in  the  United  States  carried  out 
firing  tests  of  caliber  .30  tubular  projectiles  in  1894.  A  description  of 
those  test  firings  and  an  updating  of  the  results,  performed  by  I.E.  Segal,  is 
included  in  Ref.  2.  In  the  1894  tests,  a  vulcanized  fiber  sabot  was  used  to 
seal  the  projectile  while  in  the  barrel. 

The  purpose  of  the  1894  test  was  to  compare  trajectories  of  standard  and 
tubular  projectiles.  This  was  done  by  comparing  the  vertical  drop  on  target  at 
a  given  range.  Segal's  report  indicates  that  the  drag  coefficients  computed 
from  the  1894  results  agree  closely  with  those  of  Charters  and  Thomas  [2], 
even  though  the  latter  wnrk  was  reported  fifty  years  later  (1944)  and  utilized 
more  sophisticated  experimental  apparatus. 


lumbers  appearing  in  brackets  refer  to  the  list  of  references. 


RECENT  DEVELOPMENTAL  STUDIES 

The  result  obtained  by  independent  researchers  as  to  the  relative  merits 
of  the  tubular  projectile  as  an  alternative  to  the  conventional  round  are  con¬ 
tradictory.  It  seems  as  if  there  was  as  much  research  devoted  to  discrediting 
the  tubular  projectile  as  there  was  to  investigating  its  merits.  Frank  and 
McLaughlin  [3]  have  accumulated  a  great  deal  of  data  from  various  sources  in  an 
attempt  to  "objectively  compare"  the  merits  of  the  tubular  and  "conventional 
shapes".  The  authors  concluded  that  the  tubular  has  no  particular  advantages 
over  "well  designed"  conventional  projectile  shapes.  Their  findings  disagree 
with  many  of  those  discussed  in  the  following  sections. 

Range  Tests 

Winchenbach,  Daniel  and  Edgar  [4]  conducted  range  tests  of  six  configurations 
of  tubular  projectiles  and  concluded  that  the  drag  coefficients  were  significant¬ 
ly  lower  than  the  standard  High  Explosive  Incendiary  (HEI)  projectile  of  the  same 
caliber.  Only  projectiles  of  the  same  bore  size  were  compared.  The  experimental 
models  were  constructed  from  standard  20-mm  ammunition  by  boring  various  hole 
sizes  to  obtain  "...area  ratios  (A^/A-j)  of  0.7,  0.8,  0.9,  and  1.0."  Where  At  is 
the  throat  area  and  A^  is  the  inlet  area  of  the  projectile.  The  drag  coefficients 
of  the  projectiles  with  area  ratios  of  0.8,  0.9,  and  1.0  were  less  than  half  that 
of  the  projectile  with  area  ratio  of  0.7.  The  higher  drag  coefficient  was  attri¬ 
buted  to  the  underexpanded  flow  at  the  exit  resulting  in  high  base  drag. 

Range  tests  conducted  for  the  concept  evaluation  of  the  20-mm  tubular  pro¬ 
jectile  for  the  Vulcan  Weapons  System  [5]  showed  that  the  tubular  projectile  per¬ 
formed  better  than  the  standard  HEI  round,  with  lower  drag,  shorter  time  of  flight 
for  a  given  range  (30  percent  shorter  at  1000  meters,  40  percent  shorter  at  2000 
meters)  and  similar  dispersion  characteristics  on  target. 

Recent  tests  at  the  Naval  Weapons  Center,  China  Lake,  California  of  the  Sal  1  - 
Obturated  Tubular  Projectile  (80T)  have  shown  a  definite  advantage  in  that  the  low 


drag  of  the  tubular  projectile  means  slower  retardation  of  velocity  and  hence 
a  higher  terminal  kinetic  energy  [6].  (The  BOT  is  the  object  of  this  investi¬ 
gation  and  will  be  described  in  detail  in  later  sections  of  this  report.) 

Target  Impact  Tests 

The  higher  the  energy  delivered  to  the  target,  the  greater  the  damage. 

Since  kinetic  energy  is  directly  proportional  to  the  first  power  of  mass  and 
to  the  square  of  the  velocity,  if  a  tubular  projectile  has  relatively  less  mass 
this  maybe  compensated  by  a  higher  muzzle  velocity.  However,  a  reduction  in 
mass  to  obtain  higher  muzzle  velocities  may  not  be  necessary,  depending  upon 
the  imposed  performance  criteria. 

Target  impact  studies  have  been  carried  out  with  various  projectile  con¬ 
figurations  and  target  types.  Rhethorst,  et.al.  [7]  conducted  impact  studies 
of  7.62-mm  tubular  projectiles  on  helmets.  Tests  showed  that  even  with  the 
same  energy  of  impact,  the  tubular  projectile  penetrated  further. 

Kitchen  and  Keeser  [8]  conducted  studies  for  the  Air  Force  on  the  impact 
effectiveness  of  tubular  projectiles  on  simulated  aircraft  fuel  cells.  These 
tests  were  conducted  with  steel  and  depleted  uranium  (DU)  tubular  projectiles 
and  standard  20-mm  HE^  projectiles.  The  projectiles  were  fired  at  double 
panels  at  varying  degrees  of  obliquity.  Of  the  forty-seven  tests  of  the  steel 
tubular  projectiles  which  impacted  the  target,  twenty  breached  the  rear  panel 
when  fired  at  angles  up  to  70°.  The  DU  projectiles  breached  the  rear  panel 
even  at  angles  of  85°  and  fires  were  started  in  three  of  the  five  DU  tests. 

The  standard  M56  HEI  projectiles  failed  to  breach  the  rear  panel  in  each  of 
five  firings  even  though  severe  damage  to  the  front  panel  and  two  fires  occurred. 

Brunsvold  and  Kalivretenos  [9]  conducted  a  program  to  test  the  effective¬ 
ness  of  the  tubular  projectile  against  a  simulated  cruise  missile  warhead. 

Tests  were  conducted  with  20  and  40-mm  rounds.  Results  of  those  tests  are 
classified,  and  therefore  not  presented  here.  It  was  noted  [9]  that  due  to  the 
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improved  stability  characteristics  of  the  tubular  projectile,  only  about 
half  the  spin  rate  need  be  imparted  by  the  rifling  in  the  barrel  to  obtain 
the  same  stability  as  a  standard  projectile,  with  obvious  implications  with 
respect  to  wear. 

Weapons  System  Compatibility 

For  the  tubular  projectile  to  be  practical,  it  must  be  able  to  interface 
with  existing  weapons  systems.  Reference  5  is  a  report  of  the  1978  tests 
conducted  by  the  Army  to  evaluate  the  tubular  projectile  in  the  Vulcan  weapons 
system. 

Because  the  nose  of  the  tubular  projectile  is  flat  compared  to  the  stan¬ 
dard  projectile  (See  Fig.  1),  it  "...did  not  lend  itself  to  chambering  in  the 
weapon."  Personnel  safety  required  that  the  weapon  be  remotely  operated  and 
therefore  only  the  surface-to-surface  performance  comparison  was  made. 

The  dispersion  on  target  improved  for  the  tubular  projectile  at  ranges  of 
from  700  to  2000  meters.  The  dispersion  varied  on  the  average  by  only  0.2 
mil  over  the  entire  range  spectrum  between  the  tubular  and  conventional  HE  I 
projectile.  However,  at  2000  meters,  the  tubular  was  better  by  0.4  mil.  The 
tubular  projectile  had  a  30  percent  shorter  flight  time  at  1000  meters  and  a 
40  percent  shorter  flight  time  at  2000  meters  than  the  conventional  HEI 
projectile. 

Results  of  firing  of  25 -mm  and  30-mm  tubular  projectiles  from  the  Oerlikin 
KBA  and  GAU-8/A  respectively  were  reported  in  part  in  Ref.  3  Figures  comparing 
Line-of-Sight  Penetration  Capability  vs.  Range  for  a  conventional  spinner  con¬ 
figuration  based  on  the  AR-2  shape,  a  finner  configuration  based  on  the  F10 
shape,  and  a  tubular  projectile  were  shown.  The  tubular  projectile  was  shown 
to  be  inferior  in  both  instances.  The  report  does  not  indicate  the  degree  of 
compatibility  the  tubular  configuration  has  with  the  guns  used. 

Internal  Shock  Wave  Considerations 

From  photographs  of  tubular  projectiles  in  flight  [4]  and  during  wind 
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tunnel  testing  [7],  various  shock  patterns  have  been  observed  at  the  inlet 
region  and  in  the  wake.  The  photographs  of  Figs.  2a-c  show  the  BOT  and  the 
conventional  projectiles  in  flight.  The  sabots  used  in  Refs.  4  and  7  were 
of  the  pusher  type  consisting  of  a  simple  base  plug.  The  detached  bow 
shock  seen  in  the  photograph  of  the  standard  projectile  (Fig.  2a)  is  similar 
to  that  of  the  tubular  projectile  shown  in  Fig.  2b  with  the  tubular  passage 
closed  (and  in  Ref.  4  with  the  sabot  attached).  With  no  blockage  of  the 
internal  passage  (Fig.  2c)  there  was  no  detached  bow  shock,  with  only 
attached  oblique  shocks  emanating  from  the  lip  of  the  projectile.  This  latter 
case  allows  the  possibility  of  four  other  internal  flow/shock  configurations. 

(1)  A  normal  shock  standing  at  the  entrance  to  the  projectile; 

(2)  A  normal  shock  standing  at  some  intermediate  position  in  the 
channel  ; 

(3)  A  system  of  oblique  shocks  present  in  the  channel; 

(4)  The  channel  is  devoid  of  . shocks  and  the  flow  is  shock-free 
through  the  channel . 

Reference  10  contains  a  discussion  of  the  shock  patterns  and  the 
various  conditions  leading  to  their  existence. 

PURPOSE  OF  STUDY 

The  results  of  tests  and  experiments  previously  mentioned  show,  for 
the  most  part,  definite  advantages  accruing  to  the  tubular  projectile  over 
the  conventional  projectile  primarily  in  the  area  of  lower  drag  (approximately 
1/2  to  1/3  of  conventional  [8])  which  allows  a  flatter  trajectory  and  shorter 
time  of  flight  for  a  given  range  [5].  With  reference  to  air  combat  and 
anti-missile  defense  both  from  a  Surface- to-Air  and  Air-to-Air  standpoint, 
these  features  are  definitely  worth  further  study. 

The  problem  seen  from  the  air  platform  is  that  of  tie  discarding  sabot 
or  pusher.  The  solid  [4,  5,  7,]  or  split  [8]  disk  is  effective  and  reliable. 
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30T  with  passage  blocked  resulting  in 
detached  bow  shock.  (Photograph 
courtesy  of  NWC,  China  Lake) . 


BOT  with  no  blockage  in  passage;  only  attached 
oblique  shocks  emanating  from  the  lip  of  the 
projectile.  (Photograph  courtesy  of  NWC,  China 


however,  the  possibility  of  ingestion  in  the  aircraft  engines  prohibits 
their  use.  Rhethorst  et  al . ,  in  Ref.  7  illustrates  a  number  of  sabot/ 
obturator  designs  for  use  with  the  tubular  projectile,  but  these  are  also 
of  the  discarding  type.  They  also  investigated  the  possibility  of  using  a 
consumable  sabot  [11]  which  would  be  burned  up  as  the  projectile  was  tran¬ 
siting  the  barrel  and  would  be  completely  consumed  as  it  exited  the  barrel. 

The  most  recent  development  has  been  the  BOT  which  was  designed  at 
NWC  China  Lake,  California  [6].  The  ball  obturator  has  been  bored  with  a 
hole  the  same  diameter  as  the  hole  through  the  projectile.  When  loaded, 
the  ball  is  supported  inside  the  projectile  such  that  the  axis  through  the 
hole  in  the  ball  makes  an  angle  of  approximately  90°  with  the  axis  of  the 
hole  through  the  projectile  (see  Fig.  3).  Gas  pressure  from  the  burning 
propellant  holds  the  ball  obturator  fixed  with  respect  to  the  projectile 
(as  well  as  to  the  ball).  Upon  exiting  the  barrel,  the  propellant  gas 
pressure  is  released  and,  in  its  place,  forces  due  to  aerodynamic  effects 
come  into  play.  At  that  time  a  complex  force  distribution  acts  on  the  ball 
to  cause  it  to  change  orientation  inside  the  projectile  in  such  a  way  as  to 
align  the  holes.  Inertial  forces  then  dominate  to  maintain  this  alignment 
so  that  the  projectile  remains  fully  tubular.  There  are  no  scnarating 
parts  and  the  opening  process  is  automatic. 

The  purpose  of  the  study  reported  herein  has  been  to  analytically  predict 
the  motion  of  the  ball  obturator  as  a  function  of  time.  In  support  of  this 
purpose,  an  experimental  apparatus  was  designed  and  built  to  simulate  the 
motion  of  the  obturator  in  order  to  provide  physical  evidence  of  the  accuracy 
of  the  model  and  to  gain  insights  into  the  nature  of  the  motion. 
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Figure  3 


Sketch  showing  dimensions  and  positioning  of  the  ball 
obturator  within  the  spherical  cavity. 
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ANALYSIS 


DEVELOPMENT  OF  THE  EQUATIONS  OF  MOTION 

The  ball  obturator  is  a  rigid  body  of  revolution  with  a  system  of 
coordinate  axes  fixed  to  the  ball  (body-fixed  axes)  having  its  origin  at 
the  mass  center  of  the  ball  (Fig.  4).  These  axes  are  designated  x,  y,  and 
z,  where  z  is  the  axis  through  the  hole  in  the  ball. 

To  describe  the  motion  of  the  ball  relative  to  the  projectile  an 
inertial  frame  of  reference  is  defined  with  axes  designated  X,  Y  and  Z 
with  its  origin  also  at  the  mass  center  of  the  ball.  For  this  analysis 
it  is  assumed  that  the  projectile  is  constrained  to  rotate  about  the  Z-axis. 
Because  the  motion  of  the  ball  with  respect  to  the  projectile  is  of  special 
interest  it  should  be  noted  that  components  of  the  ball  angular  velocity 
with  respect  to  the  X  and  Y  axes  are  also  relative  to  the  projectile, 
whereas  the  component  with  respect  to  the  Z-axis  must  account  for  the  pro¬ 
jectile  spin. 

The  motion  may  be  described  by  Euler's  Modified  Equations  of  Motion 
for  a  rigid  body  of  revolution  about  a  fixed  point.  These  are  given  in 
Thomson  [12]  as: 

Adi  +  (C-A)u  &)  =  M  (la) 

1  2  3  1 

Aw  +  (A-C)ci)  to  =  M  (lb) 

2  13  2 

Cw  =  M  (lc) 

3  3 

where 

A  =  Mass  moments  of  inertia  about  principal  axes 

perpendicular  to  the  z-axis. 

C  =  Mass  moment  of  inertia  about  the  principal 

axis  through  the  hole  (z-axis). 


Figure  4.  Illustration  of  coordinate  systems  and  symbols  used  in  the 
analysis. 
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ci)  ,  go  ,  co  a  Angular  accelerations  about  the  body-fixed  axes 

12  3 

x,  y,  z  respectively. 

a)  ,  a)  ,  a)  =  Angular  velocities  about  the  body-fixed  axes. 

12  3 

M  ,  M  ,  M  *  External  moments  acting  about  the  body-fixed 

12  3 


The  mass  moments  of  inertia  for  the  particular  case  at  hand 
(sphere  of  radius  R  with  concentric  hole  of  radius  r)  are  given  by: 


47rpRs 

=  ~T5 
_  4icR5 


[^(1-r2) 2  (3r2  +  2)  +  (1  -  r2)s'2 


C  *  (1-r2)3'2  (3?2  +  2)  (3) 

where  r  =  r/R. 

The  symbol  X  is  introduced  here  to  denote  the  relative  difference 

between  the  two  moments  of  inertia.  That  is: 

,  .  C-A 

x  ‘  ~r 

From  Eqs.  (2)  and  (3)  it  is  seen  that 

x  =  £r  (4> 

and  since  0  -  r  -  1,  X  lies  in  the  range  0  -  X  -  1. 

With  X  so  defined,  Eqs.  (1)  may  be  rewritten: 


CD  +  XCD  CD 

=  H 

/A 

1  2  3 

i 

• 

CD  -  XCD  CD 

=  M 

,/A 

2  1  3 

2 

• 

CD 

=  M 

/c 

3  3 

The  orientation  of  the  ball  at  any  time  may  be  described  by 
three  angles;  9,  $  and  i|».  These  are  known  as  Euler  angles  and  their 
relationships  to  the  coordinate  axes  already  described  may  be  seen  in 
Fig.  4.  The  components  of  angular  velocity  coj ,  cd  and  cd  in  terms  of 
the  Euler  angles  are  given  as  [12]: 
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Similarly,  if  the  multiplication  is  reversed  and  the  results  are 
subtracted,  we  have 

n  “  ”  “ 

9  =  (R  /A)cos$  -  (M  /A) s i nd>  -  ^sin8  [(X  +  1)  $  +  Xijicose] 

1  2 

Simply  rewriting  Eq.  (7c)  we  have 

t*  -  n 

0  =  R  /C  +  ijesine  -  i|>cos9 

3 

These  are  the  basic  governing  equations  for  the  motion  of  the  ball 
in  terms  of  the  Euler  angles.  They  are  second  order  ordinary 
differential  equations,  nonlinear,  and  strongly  coupled.  Before 
proceeding  further,  we  shall  introduce  the  convenience  of  non- 
dimensional  ization  of  the  time-dependent  variable.  This  is  ac¬ 
complished  through  the  use  of  ojpo,  the  initial  rate  of  spin  of  the 
projectile.  The  equations  of  motion  are  unchanged  in  form  except 
for  the  terms  involving  the  external  moments. 

Thus,  with  the  definitions 

M  =  R  /Au  2,  M  =  M  /Aco2  M  =  fl  /Co  2 
i  i'  po  2  2'  po  3  3  po 

we  have 

\psine  -  Misin4>  +  M2cos<J>  +9  [(X  +  1)  4>  +  (X-l)$cos6]  (8) 

0  =  MiCOS<fi  -  M2sin<f>  -  $sin0  [(X  +  1)<J  +  Xi^coss]  (9) 

>  =  M3  +  tji9sin8  -  ij/cos0  (10) 

—  ••  n 

where  ij  =  ^/a^0,  rp  =  4>/upo2 ,  etc. 

The  integration  of  Eqs.  (8)  -  (10)  proceeds  from  a  condition 
in  which  the  ball  is  fixed  to  the  projectile  and  spinning  with  it. 
The  initial  values  of  the  precession  angle  \p  and  the  spin  angle  <j> 
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are  arbitrary  and  are  set  to  zero  for  convenience.  The  appropriate 


initial  conditions 

are,  therefore: 

O 

CD 

II 

O 

CD 

9(0)  =  0 

(U) 

4»(0)  «  0 

$(0)  =  0 

O 

II 

o 

$(0)  =  1 

The  central  problem  yet  to  be  addressed  is  that  of  the  specification 
of  the  external  moments  ft  i  ft  ,  and  M  or,  in  non-dimensional  form, 

12  3’ 

M  ,  M  ,  and  M  . 

1  2  3 

EXPRESSION  OF  THE  EXTERNAL  MOMENTS 

The  interaction  of  the  ball  with  the  cavity  in  the  spinning 
projectile  is  dominated  by  torques  arising  from  the  relative  motion 
between  the  two  bodies.  These  torques  will  be  opposite  to  the  direction 
of  the  relative  velocity  so  that  we  may  begin  by  specifying  that: 

M  =  -1  M  (12) 

where  lr  =  <or/|(ur|  *  unit  vector  in  the  direction  of  the  relative 

a/  *+* 

angular  velocity 

We  have  considered  torques  due  to  both  fluid  shear  and  sliding  friction. 


Nakabayashi  [13]  has  recently  reported  the  results  of  a  series  of 
experiments  to  determine  the  fluid  shear  acting  in  the  gap  between  two 
concentric  rotating  spheres.  Although  the  details  of  the  flow  are 
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extremely  complex  [14],  the  results  of  Nakabayashi  show  that  the 
data  are  well  correlated  by: 


Cm  =  87r(l+  e)/(3eRe)  for  laminar  flow  (13) 

and  Cm  =  0 . 053 [  1  +  (7/4)e]/(e[l  -  (3/2)e]Re°  ‘ 25}  for  turbulent  flow, 
where  Cm  =  M^/(pR5(Jr2) ,  e  =  (h/R)<<l,  and  the  Reynolds  number  is  given 
as: 


Re  =  R2cur/v. 

From  [13]  it  is  also  possible  to  estimate  a  transition  Reynolds 
number: 

Ret  =  70e-1-6 

In  our  experiments  (both  real  and  numerical)  we  have  found  that  relative 
angular  velocities  in  the  range  of  interest  are  well  below  those 
necessary  to  stimulate  turbulent  flow  in  air  so  that  Eq.  (13)  is 
adequate  for  our  purposes.  When  Eq.  (13)  is  cast  in  a  suitable  form  we 
have: 


Mf  =  8ir(l  +  e)uR3ir/(3e)  (14) 

Torque  Due  to  Sliding  Friction, 

In  the  case  of  sliding  friction  the  magnitude  of  the  applied  torque 
may  be  written 


»s  '  “s  Fn  rm 


(15) 


where  Fn  is  the  normal  component  of  the  resultant  of  all  forces  acting 
upon  the  ball  and  rm  is  the  moment  arm  extending  from  the  point  of  con¬ 
tact  and  directed  normal  to  the  vector  of  relative  angular  velocity. 

Since  the  point  of  contact  must  be  on  the  surface  of  the  ball, 
the  position  vector  of  this  point,  P,  must  have  a  magnitude  of  the 
ball  radius  R.  Figure  5  illustrates  the  geometry  for  cases  to  be  dis¬ 
cussed  in  which  multiple  points  of  contact  are  distributed  along  a  line. 

If  y  and  B  are  the  location  angles  of  a  point  of  contact,  as  illustrated 
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in  Fig.  5,  then  the  position  vector  nay  be  written 
P  =  R(e  sinBcosy  +  e  sinSsirry  +  e  co-:J) 

-V  ■vx  \z 

where  e  ,  e  ,  and  e  are  unit  vectors  for  the  inertial  coordinates. 
%x  n/  ^ z 

The  vector  of  relative  angular  velocity  may  be  written, 
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(19) 


and  the  magnitude  of  the  moment  is 

Msg  =  MSWR[l-(<urx/wr)  2]lj 

Which  is  applied  opposite  to  the  relative  velocity  vector  in  accordance 
with  Eq.  (12). 

Axial  aerodynamic  load.  The  aerodynamic  load  acting  upon  the 
ball  will  be  a  complex  function  of  the  ball  orientation.  During  the 
opening  instants,  however,  when  the  ball  blocks  the  tubular  projectile, 
the  resultant  of  the  pressure  forces  on  the  ball  will  act  along  the  axis 
of  the  projectile.  This  interval  is  of  importance,  of  course,  since  it 
will  constitute  a  significant  part  of  the  ball-opening  delay.  At  lower 
values  of  9,  when  there  is  significant  flow  through  the  ball,  the  specif¬ 
ication  of  the  net  aerodynamic  force  will  require  further  detailed  analysis. 
It  is  important  to  note  further  that  when  the  ball  is  partially  open  there 
may  be  a  significant  moment  acting  upon  the  ball,  as  well  as  a  centrally 
directed  force.  An  understanding  of  the  nature  of  the  aerodynamic  force 
under  partial ly-open  conditions  is  a  major  goal  of  continuing  research. 

For  the  case  of  an  axially  oriented  applied  force,  the  ball  will 
be  forced  against  the  aft  rim  of  the  hole  in  the  projectile.  The  exact 
number  and  location  of  the  contact  points  is  difficult  to  predict  and 
will,  in  fact,  depend  upon  projectile  manufacturing  variations.  However, 
any  effect  of  arbitrariness  in  contact  point  specification  is  "blurred" 
by  the  high-speed  rotary  motion  and  errors  resulting  should  be  random 
and  therefore  of  little  influence  on  the  overall  motions.  This  hypo¬ 
thesis  has  been  verified  by  computer  experiments.  In  selecting  the  number 
of  contact  points  and  their  location,  we  have  assumed  the  first  stable 
position,  that  is,  three  contact  points.  In  addition,  we  have  assumed 
that  these  points  are  symmetrically  dispersed  along  the  line  of  contact. 


The  axial  aerodynamic  force,  F  .  is  evenly  balanced  among  the 

a 

contact  points  so  that  at  each  point  the  normal  force  is  given  by 

F  =  ^a  cosS 
n  3 

and,  since  the  contact  points  lie  along  the  rim  of  the  projectile 
hole,  8  =  sin'1  r/R.  (see  Fig.  5).  The  moment  due  to  sliding  friction 
in  this  case  is  given  by 

Bsa  *  Ms(fa/3)co56  Mr,).  (20) 

where  8  =  sin'1 (r/R)  and  the  three  values  of  rm  are  given  by  Eq.  (18) 
with  the  angle  y  set  at  0,  2tt/3,  and  4tt/3. 

The  aerodynamic  force  most  appropriate  for  use  in  Eq.  (20)  will 
be  that  due  to  a  pressure  distribution  over  the  ball  and  behind  a  normal 
shock  standing  at  the  projectile  entrance,  (see  sketch  below) 


The  net  force  on  the  ball  will  be 

F  =  7rr2  ( P  -  P  ) 
a  2i 

Since  the  actual  base  pressure  Px  is  much  less  than  P2  and,  in  any  case, 
is  unknown,  we  assume  ambient  pressure  here.  That  is,  Pj  *  PM.  The  face 


P2  will  lie  somewhere  between  the  static  and  the  stagnation  pressure 
existing  behind  the  shock.  The  difference  is  small  (about  7 *„)  in  the 
range  of  Mach  numbers  of  interest  (  =  3  to  5)  and  will  likely 

approach  the  static  pressure  because  of  spillage  behind  the  shock  and 
recirculation  within  the  cavity  ahead  of  the  ball;  that  is,  there  will 
be  little  or  no  pressure  recovery  in  the  cavity.  The  assumption  of 
static  pressure  is  adopted  for  these  reasons  as  well  as  the  fact  that 
this  is  a  conservative  assumption  in  terms  of  ball  opening  times  (higher 
pressures  lead  to  shorter  opening  times).  Thus,  the  approximate  ex¬ 
pression  for  the  pressure  P2  is 

p2  =  p„  ti  +  Hr  (mj  -  i)] 

and  for  the  aerodynamic  force, 

Fa  =  Trr^tj^M^2  -1)]  (21) 

Equations  (20)  and  (21)  provide  the  necessary  information  for  determining 
the  torque  due  to  an  axial  aerodynamic  load  and  this  torque  is,  again, 
applied  in  a  direction  opposite  to  that  of  the  relative  velocity  vector. 
Consideration  of  the  Relative  Angular  Velocity  \fector 

The  three  body-fixed  components  of  the  relative  velocity  vector  are 
given  by  see  Eqs.  (5)  : 

a)jr  =  (4»  -  a)p)sin0sin<j>  +  0cos<f> 

&>2r  =  (4>  -  (Up)sin9cos (j>  -  0sin4>  (22) 

“jr  =  (*!>  *  cop)cose  +  $ 

According  to  Eq.  (12),  then,  the  components  of  the  external  moments 


may  be  written: 


°po  “r 


upo  r 

R 


[(4»  -  oj  ) s i n6s i n4>  +  0cos<b] 


[(<£  -  6jp)sinecos<}>  -  8sin<j>] 


c^7Sr  [(*  -  “p)cos6  * 
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The  terms  involving  these  moments  in  the  equations  of  motion  may 


now  be  written 

M  sin<£  +  M  cos^  = 
1  2 

M  cos<J>  -  M  sin*  = 
1  2 


A“po2“r 

-M 

A“po*®r 


I(4>  -  cop ) s i ne ] 

0 


=  as  in  Eq.  (23) . 


(24) 


The  value  of  R  in  these  expressions  are  given  according  to  the  table 
below: 

Source  of  moment  Applicable  Equations 

(14) 


Fluid  Shear, 

Sliding  Friction 
Gravity  only,  Msg 
Axial  Aerodynamic,  R 


sa 


(19) 

(20,  21) 


It  is  appropriate  to  combine  the  fluid  shear  moment  with  either 
of  the  sliding  friction  moments.  In  most  cases  analyzed,  however, 
fluid  shear  moments  are  several  orders  of  magnitude  smaller  than 
those  due  to  sliding  friction.  Since  the  gravity  and  axial  aerodynamic 
loads  are  each  dependent  upon  specified  points  of  contact,  it  is  not 
appropriate  to  combine  them  without  determining  a  point  of  contact  that 
is  correct  for  the  combination.  In  cases  of  practical  interest,  however, 
the  aerodynamic  load  is  either  absent  or,  if  it  is  present,  it  completely 
dominates  any  effects  due  to  gravity.  Consider,  for  instance,  a  sea  level 
launch  under  the  following  conditions: 
r  *  4.7  nm 
R  =  7.9  run 
K,  =  3.0 

Under  the  approximation  that  the  moment  arms  are  about  equal  to  the  sphere 


radius  for  both  gravity  and  aerodynamic  loads,  the  ratio  of  moments 
for  the  two  cases  is 


Msg  =  W 

"Sa  fa 


For  the  stated  conditions  with  a  steel  ball,  W  =.08  N  (about  0.4  lb.) 
and  F  =  65  N  (about  15  lb.)  so  that  the  aerodynamic  effect  is  clearly 

a 

dominant. 

It  should  be  noted,  in  closing  this  section,  that  the  main 
weakness  in  the  analysis  is  due  to  uncertainties  in  specifying  the 
aerodynamic  load.  If  this  load  (force  vector  or  centrally  applied 
force  plus  moment)  is  known,  the  methods  given  to  predict  its  affect 
upon  the  motion  are  quite  general. 


EXACT  SOLUTIONS  FOR  TYPICAL  SYSTEM  PARAMETERS 


Equations  (8)  -  (10),  subject  to  the  initial  conditions 
[Eq.  (11)]  have  been  solved  for  both  the  gravity-only  case  and  the 
case  of  an  axial  aerodynamic  load.  The  gravity-only  case  is  of  special 
interest  because  it  provides  theoretical  predictions  for  the  experimen¬ 
tal  program  (described  in  a  subsequent  section).  The  aerodynamic  load 
is,  of  course,  the  practical  case  for  the  tubular  projectile. 

The  appropriate  equations  were  coded  for  solution  on  the  digital 
computer.  The  general  flow  of  the  solution  method  is  as  follows: 

1.  Set  initial  conditions  and  specify  system  parameters. 

2.  Calculate  necessary  constants  (ball  weight,  moments  of  inertia, 
flight  level  pressure,  etc.). 

3.  Calculate  a)  ,  co  ,  <o  from  current  Euler  angles  [Eqs.  (5)]. 

i  %  3 

4.  Calculate  body-fixed  components  of  relative  angular  velocity 
[Eqs.  (22)]. 

5.  Convert  relative  angular  velocity  components  to  inertial  axes 

(ai  ,  a>  ,oj  ) . 

\  rx»  ry  rz; 

6.  Determine  moment  arms  and  external  moments  as  required  (see 
table,  p.  22). 

7.  Determine  moment  terms,  for  use  in  equation  of  motion  [Eqs.  (24)]. 

8.  Determine  ip,  9,  and  <p  [Eqs.  (8)  -  (10)]. 

9.  Integrate  for  new  values  of  4>,  0,  and  4s  0,  and  <j>. 

10.  Return  to  step  3,  until  0  is  sufficiently  close  to  zero. 

The  method  of  integration  employed  (step  9  above)  was  a  fourth-order 
Runge-Kutta  scheme  with  a  step-size  automatically  adjusted  to  control 
numerical  precision. 
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.  nsideration  of  Projectile  Angular  Velocity 
Two  system  constraints  have  been  investigated  in  the  analysis. 

In  the  first  case  we  have  invoked  the  constraint  of  constant  pro¬ 
jectile  spin  rate  and,  in  the  second  case,  we  have  required  a  con¬ 
stant  system  (ball  plus  projectile)  kinetic  energy.  Since  the  in¬ 
ertia  of  the  ball  is  small  relative  to  that  of  the  projectile  (a 
factor  of  about  20)  the  decrease  in  projectile  spin  rate  necessary 
to  maintain  constant  system  kinetic  energy  during  the  ball  motion 
is  rather  small  (typically  about  2%)  and  the  two  cases  are  not  sig¬ 
nificantly  different.  The  computer  solutions  shown  here  (so-called 
"exact"  solutions)  are  for  the  case  of  constant  system  kinetic  energy. 
Gravity  Load. 

The  system  parameters  for  this  case  were  set  to  correspond  to 
the  conditions  of  the  experiments  to  be  described  later.  These  para¬ 
meters  are  summarized  in  Table  I. 

Figures  6  and  7  illustrate  the  nature  of  the  motion  predicted  by 
the  exact  solutions  for  a  representative  set  of  the  physical  parameters. 

In  Fig.  6  the  angle  of  tilt  (9)  relative  to  the  cavity  spin  axis  is 
shown  as  a  function  of  time  for  various  rates  of  spin.  At  relatively 
high  rates  of  spin  the  motion  appears  to  be  quite  linear  with  a  nearly 
constant  rate  of  nutation  (9)  developed  early  in  the  motion.  At  lower 
rates  of  spin,  however,  there  is  an  initial  period  during  which  the  mo¬ 
tion  is  nearly  stable,  followed  by  a  rapid  decrease  in  the  tilt  angle. 

The  time  required  to  reach  a  given  value  of  9  is  therefore  a  function  of 
the  angle  itself  and  is  illustrative  of  the  complex  interactions  between 
the  inertial  tendencies  of  the  ball  and  the  dissipative  torques  due  to 
friction.  Figure  7  is  included  to  show  the  variations  in  the  Euler  angle 
rates  for  a  typical  case.  When  the  ball  is  released  its  rate  of  precession 
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TABLE  I 


Ball  Characteristics  and  Experimental  Conditions  -  Gravity  Load 


Ball  Designation 

Mod  0 

Mod  I 

Physical 

Characteristics 


R,  m 

7.899  x  10"3 

7.925  x  10'3 

r,  m 

4.722  x  10‘3 

1.984  x  10‘3 

C,  Nms2 

3.126  x  10'7 

4.029  x  10"7 

A 

0.410 

0.077 

P,  kg/m3 

7.67  x  103 

7.75  x  103 

Experimental 

Conditions 

u,  Ns/m2 
Ps  (note  l) 
e 

90»  deg 
9*,  deg 
R,  Nm  (note  2) 


1.917  x  10‘5 
0.51 

6.43  x  10’5 

87.25 

36.71 

2.08  x  10'4 


1.917  x  10"5 
0.30 

6.43  x  10"5 

89.5 

14.50 

2.18  x  10-4 


^Empirically  determined,  see  text. 
2 

Computed  from  Eq.  (39). 


» 
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Figure  7.  Euler  angle  rates  (non-dimensional)  for  typical 
case.  Gravity  load  only,  ball  Mod  0,  <*>  =  400 


(4»)  increases  above  its  initial  non-dimensional  value  of  unity  and, 
concurrently,  a  retrograde  spin  is  developed.  For  the  conditions  of 
Fig.  7  (<Up0=  400s*1 ),  the  mean  value  of  nutation  rate  is  fairly  constant 
with  acceleration  rates  that  are  initially  large  and  rapidly  damped. 
Additional  results  for  case  of  gravity  load  are  discussed  in  con- 
juntion  with  the  linear  approximation  and  experimental  results. 

Axial  Aerodynamic  Load. 

Figures  8  and  9  present  similar  results  for  the  motion  under  an 
aerodynamic  load.  The  system  parameters  for  these  figures  are  those 
of  the  "baseline  case"  listed  in  Table  II. 

The  nature  of  the  response  is  seen  to  be  similar  to  that  of  the 
gravity-only  case.  In  particular,  the  existence  of  a  hovering  mode  is 
apparent  at  low  projectile  spin  rates.  This  phenomenon  points  up  some 
design  limits  and  is  the  subject  of  discussion  in  a  subsequent  section. 
Sensitivity  studies  were  conducted  for  the  aerodynamic  load  and  these 
also  are  discussed  later. 
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TABLE  II 

Ranges  of  Design  Variables  -  Aerodynamic  Load 


Design  Variable 

Minimum 

Basel ine 

Maximum 

Figure  No.1 

R,  mm 

base! ine 

7.90 

39.5 

24 

A,  Nms2 

2.24xl0"7 

7.00x10"* 

r/R 

0.400  7 

0.598 

0.800 

25 

A,  Nms2 

3.20x10' 

2.24xl0"7 

l.OOxlO"7 

\ 

0.19 

0.41 

0.69 

p,  kg/m3 

3.10xl03 

7.75x10s 

1.24x10* 

26 

A,  Nms2 

8.96x10"® 

2.24xl0"7 

3.58X10"7 

^s 

0.71 

0.855 

1.000 

27 

M 

00 

3.00 

3.50 

5.00 

28 

P  .  atm 

CO 

0.1 

1.0 

basel ine 

29 

a..,  s"1 
po 

6.91x10s 

1.20x10* 

1.36x10* 

8 

9o>  deg 

70 

85 

baseline 

8,  10 

1  Figures  cited  are 
investigated  over 

those  in  which  the  given 
a  range  of  values. 

design  variable 

i  s 

i 

* 


I 
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LINEAR  APPROXIMATIONS 


From  Figs.  8a  and  8b  it  will  be  noted  that  there  is  a  marked  change 
in  the  nature  of  the  ball  motion  at  a  projectile  spin  rate  of  about 
7300  s"1  (70,000  rpm).  At  rates  above  this  value  (Fig.  8a)  the  ball  motion 
appears  to  be  quite  linear;  that  is,  the  response  of  the  ball  to  the  friction 
al  torques  approximates  one  of  a  second-order  linear  system  responding  to  a 
combined  ramp  and  step  input  (both  of  which  are  negative).  Such  a  response 
is  given  by 

9  ■  0  +  D  +  R(t  -  2?)...  (25) 

+  exp(-st)  {Rsin[(l-?2)lst  +  $R]  -  Dsin[(l -c2)^  +  $D)]}  /(l-^2)*5 
where  <1>D  =  $R/2  =  tan' 1  [(1  -s2)^] 

The  parameters  of  the  response,  to  be  derived  below,  are  given  by: 

5  =  M/2,  D  =  -(A/2)sin20o,  and  R  =  MD/U  +  1) 

The  method  leading  to  the  linear  approximation  is  a  perturbation  upon  the 
equilibrium  solution  to  Eqs.  (8)  -  (10)  under  torque-free  conditions. 

DEVELOPMENT  OF  THE  PARAMETERS  OF  THE  LINEAR  RESPONSE 

Equation  (25)  satisfies  the  initial  conditions  on  0  and  is  the  solution 
to 

9  +  2^9  +  9  =  9Q  +  D  +  Rt  (26) 

From  the  equations  of  motion  for  9  and  the  associated  relation  for  external 
moment  [Eqs.  (9)  and  (24)],  a  comparison  of  the  coefficient  of  the  tilt  ve¬ 
locity,  9  ,  leads  to  the  damping  coefficient: 

C  =  M/2  (27) 

where  M  =  M/Aui  uL . 

po  r 

With  this  definition  of  4  ,  combination  of  Eqs.  (9)  and  (26)  leads 


to 

(0  -  0  )  -  (0  +  Rt)  -  ^s1n0[(A  +  1)|  +  Ai^cos0]  (28) 

The  response  specified  in  Eq.  (25)  indicates  that  the  function  described 
in  Eq.  (28)  is  oscillatory  in  nature  with  the  mean  value  decaying  exponenti¬ 
ally  from  an  initial  value  of  -0  to  an  asymptotic  final  value  of  -2R£.  That 
is,  the  mean  value  (about  which  the  oscillations  occur)  is  given  by: 

C(9  -  0O)  -  (0  +  Rt)]mean  =  -2R?  -  (D  -  2Re)exp(-ct) 

In  the  limiting  case  of  vanishingly  small  external  moments  this  expression 
simplifies  to  a  straight-line  response  with  slope  R  so  that  an  evaluation 
of  this  rate  of  decay  should  emerge  from  an  investigation  of  the  solution  of 
the  equations  of  motion  for  the  case  of  zero  external  moments. 

If  the  external  moments  are  removed  from  Eqs.  (8)  -  (10),  a  solution 
of  the  following  form  is  obtained; 

(A  +  1)$0  +  ACOS0O  «  0,  0Q  »  0,  ^  -  1  (29) 

The  subscript  (  )  is  used  here  to  indicate  quantities  pertaining  to  this 
solution.  Perturbations  to  this  solution,  designated  by  subscript  (  ),  are 
defined  such  that 

'P  =  iP0  +  M’I'p 

<p  =  <pQ  +  M4>p  (30) 

0  =  0  +  M0_ 

O  p 

The  base  quantities  (subscript  o)  must  satisfy  Eq.  (29)  and  the  perturbed 
quantities  (no  subscript)  must  satisfy  the  initial  conditions,  Eq.  (11). 

The  appropriate  initial  conditions  for  the  perturbations  are,  therefore 

*  (0)  *  0 

♦  (0)  =  -4>0/M  =  -Acos0o/M(A  +  1)  (31) 

0p(O)  -  9p(0)  =  0 
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When  Eqs.  (29)  and  (30)  are  inserted  into  Eqs.  (8)  and  (10),  the  fol¬ 
lowing  approximate  relationships  are  obtained: 

MifpSineo  s  Ti  -  M6pCos90  +  0(M2) 

MS  =  T0  +  M§  sin0„  -  Mi|>  cose  +  0(M2) 

p  c  p  0  p  0 

where  the  torque  terms  and  are  obtained  from  Eqs.  (24),  (29),  and  (30) 
T1  =  -M(M^psin0o)  =  0(M2 ) 

T2  *  -MV(X  +  15  +  0(M2) 

Neglecting  terms  of  order  M2  and  smaller,  we  have 

Vfn9° '  •9pcos9o 

fp  -  40/(i  ♦  i)  ♦  0p«ceo 


The  first  integrals  of  this  set  yield 
=  -6  cote 

p  p  0 

Ip  3  40t/(x+  i)  +  epcsceo  -  4>q/m 

or,  from  Eqs.  (30) 

<ji  *  1  -  M6  cote 


3  -M<pQt/(X  +  1)  +  M0pCSC0o 


(32) 


(33) 


When  Eqs.  (29)  and  (30)  are  substituted  into  Eq.  (28)  we  have,  to  the  first 
order: 


1 


Msine  [(X  +  1)<b  +  XiL  cose  -  X0  sin0  ]  =  M6  -  (D  +  Rt) 

0  p  p  0  p  0  p 

and  with  Eqs.  (32)  this  becomes 

(X  +  l)sin6o[M4>ot/(X  +  1)  +  $  ]  -  D  +  Rt 

It  follows  that: 

0  =  (X  +  l)<j>osin0o  *  -(X/2)sin20Q  (34) 

R  3  MD/(X  +  1) 
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In  order  to  evaluate  the  non-dimensional  moment  parameter,  M,  it  is 
necessary  to  develop  approximate  relationships  for  and  for  the  moment 
arm  r^  to  be  used  in  the  previously  derived  expressions.  Under  the  con¬ 
straints  of  the  linearization  described  here  (M«l  and  |D|«1),  the 
relative  angular  velocity  is  approximated  by  [Eqs.  (22)]: 

u>r  ~  w3r  ~  ^o  =  '  Xcos0</(X  + 

The  principal  component  of  the  relative  angular  velocity  is  therefore 
along  the  spin  axis  of  the  ball  and,  since  the  ball  is  initially  near 
9q  =  -rr/2  (required  for  (0(«1)  the  components  along  the  inertial  axes  are, 
approximately, 

“rx^r  *  sin*o 
uryA  a  cos^o 
“rzA  =  0 

Since  i|j0  is  a  function  of  time  (even  though  ip0  *  constant  =  1  in  the 
linearization),  suitable  mean  values  for  the  moment  arm  relationships  must 
be  obtained.  In  the  case  of  a  gravity-only  load  we  have 

rmg  =  R[1  -  {aVx/a)r)2]!S=  R>C0SV 

the  absolute  value  being  necessary  to  preserve  the  positive  nature  of  the 
moment  arm.  For  the  mean  value, 

i  r2n 

r  =  7T-  I  r  dip 
mg  2tt  J  mg  yo 

0 

we  have 

i  fir  ?r  A2 

f»3  *  H  ;0  Rlcos*ol  %  ‘  irjf0  cos*o  % 

and 

?mg  =  (2/ir)R  (36) 

\ 

For  the  case  of  aerodynamic  loads,  we  have  from  Eq.  (18) 
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rma  =  '  s1‘n2ssin2(^0  +  Y )  I*5 

where  the  selected  values  of  locate  the  contact  points  and  S  =  sin_1(r/R). 


In  the  mean,  then 


?ma  ar/0  n  -  sin2Ssin2(^0  +Y)]15  d^o 


rma  *  ( 2/ir) R  E(8)  (37) 

where  E(B)  is  the  Complete  Elliptic  Integral  of  the  second  kind. 

With  these  values  for  the  moment  arms  it  remains  to  evaluate  the  line¬ 
arized  forms  of  the  external  moments.  In  the  case  of  fluid  friction  we 
have,  from  Eq.  (14) , 

Mf  *  [8tt(1  +e)yR3ir/(3b)]/Aira)po 


Mf  «  SttH  +  c  )uR3/(£Aupo)  (38) 

Note  that  no  linearization  is  required  for  this  non-dimensional  moment  para¬ 
meter.  For  the  case  of  gravity-only  loads,  from  Eq.  (19) 

Msg  -  luM(Z/n)]/fcrupo 


%  *  2»SWR(A  +  U/7rAAu)po2COS0o  (39) 

For  the  aerodynamic  load  Eq.  (37)  gives,  for  three  contact  points  symmetrically 


disposed, 


iI1?m1  «  3(2/ it) R  E(S) 


and  from  Eq.  (20) 


Msa  *  usFaR(2/7r)E(S;cos8/Ad)ru)po 


M  *  2y  F  R(\  +  1  )E(8)cos8/ttA\u)  2cos0 
sa  s  a  po  0 
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The  fluid  friction  term,  Mp  may  be  added  to  either  the  gravity  term,  M  , 
or  the  aerodynamic  term,  M.  to  obtain  M  and  the  related  parameters 

5a 

[S,  R,  and  D  from  Eqs.  (27)  and  (34)]. 

RESULTS  COMPUTED  FROM  THE  LINEAR  APPROXIMATIONS 

Figure  10  illustrates  the  linear  approximation  in  comparison  with  the 
exact  solution  for  two  initial  values  of  8  under  fluid  friction  and  gravity 
load.  (As  mentioned  previously,  the  effects  of  fluid  friction  are  generally 
negligible.  This  observation  also  applies  in  the  linearization.)  Figure  11 
gives  similar  results  of  the  case  of  an  axial  aerodynamic  load  on  the  ball. 

It  is  seen  from  both  figures  that  good  approximations  to  the  exact  solution 
are  obtained  from  the  linearization.  The  approximation  departs  from  the 
exact  solution  as  time  passes  and  the  ball  angle  decreases.  This  is  to  be 
expected  since  the  terminal  phase  of  the  ball  motion  is  not  adequately 
modeled  in  the  linearization.  Nevertheless,  the  linear  model  provides  valu¬ 
able  insights  into  the  nature  of  the  motion  and  its  most  influential  parameters. 
It  should  be  noted  again  that  at  the  lower  angles,  when  an  observer  would 
see  a  sizable  opening  through  the  projectile  when  sighting  along  its  axis, 
both  theoretical  models  are  subject  to  the  possibility  of  considerable  error. 

The  is  especially  true  in  the  aerodynamic  case  because  of  uncertainties  in 
the  load  description  at  partlally-open  conditions. 

The  results  of  the  linear  approximation  are  discussed  further  in  con¬ 
nection  with  the  experimental  program  to  be  described  in  the  next  section. 


Figure  10.  Comparison  of  approximate  solutions  with  exact  solutions 
Gravity  load  only,  ball  Mod  0,  u  =  400s"1. 
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Comparison  of  approximate  solutions  with  exact  solutions. 
Baseline  projectile,  upo  3  1 . 20x1 0U  s"‘. 
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EXPERIMENTATION 


EXPERIMENTAL  APPARATUS 

An  experimental  system  was  designed  to  simulate  the  spinning  BOT 
under  the  influence  of  a  gravity  load.  The  apparatus  consisted  of  a 
compressed-air-driven  spin-up  rig,  optical-timing  mechanism,  air  manifold 
and  associated  piping,  tubing  and  electronics. 

The  obturator  was  fabricated  from  a  standard  5/3-in.  diameter  chrome- 
steel  (52100)  bearing  ball.  The  ball  was  annealed  to  allow  machining  and 
then  bored  along  a  diametrical  axis.  The  bored  ball  was  then  mounted  in  a 
bakelite  metalographic  specimen  mount.  The  mounted  ball  was  placed  in  a 
milling  machine  and  a  flat  was  machined  in  the  bakelite  at  a  specified  angle 
relative  to  the  axis  through  the  hole  in  the  ball.  This  flat  was  then  used 
as  the  polishing  plane  for  metalographic  specimen  preparation.  A  small  flat 
spot  was  polished  on  the  obturator  to  provide  a  highly  reflective  surface  at 
a  known  orientation  relative  to  the  z-axis  of  the  obturator  (Fig.  12). 

After  removal  from  the  bakelite  mount,  the  obturator  was  placed  in  a 
three-piece,  lucite  housing  (Fig.  13).  The  mating  ends  of  the  two  hollow  in¬ 
ner  cylinders  were  each  machined  with  a  5/8-inch  spherical  end  mill  to  a  depth 
of  approximately  5/16-inch.  When  mated,  a  spherical  cavity  was  formed  to 
accomodate  the  obturator.  The  third  cylinder  was  press  fit  over  the  others 
after  the  obturator  was  inserted  to  insure  alignment  of  the  inner  pieces  and 
provide  rigidity  (Fig.  14).  The  ends  of  this  lucite  assembly  were  then  press 
fit  into  aluminum  end-pieces  similar  to  those  in  Fig.  14.  These  end-pieces 
served  as  the  shaft  for  the  bearings  and  one  also  served  as  the  prime  mover 
(bucket  wheel)  for  the  apparatus.  The  shaft  rotated  in  two  ball  bearings 
mounted  in  aluminum  pillow  blocks  aligned  on  a  rigid  pedestal.  The  prime  mover 


Obturator  and  components  of  lucite  obturator 
housing  with  mating  ends  of  inner  cylinders 
machined  to  house  the  obturator. 


Assembled  obturator  housing  with  aluminum 
end  pieces. 


was  a  bucket  wheel  machined  from  a  solid  aluminum  disk  (Fig.  15). 

The  bucket  wheel  was  driven  by  compressed  air  supplied  from  an  install¬ 
ed  system  through  an  air  filter  to  a  manifold  and  then  through  one  of  two 
Model-10  Kendall  pressure  regulators.  From  the  regulator,  the  air  passed 
through  a  flexible  tube  to  a  tee,  each  leg  of  which  supplied  a  nozzle.  These 
nozzles  were  mounted  opposite  one  another  on  the  pillow  block  in  such  a  way 
as  to  allow  the  air  jets  to  impinge  upon  the  bucket  wheel  to  cause  rotation 
(Fig.  16).  The  speed  of  rotation  was  sensed  by  a  Bentley  Nevada  Proximitor, 

Model  3100N,  which  was  mounted  above  the  bucket  wheel  to  detect  the  passage 
of  each  point  on  the  wheel.  The  pulses  generated  by  the  sensor  were  counted, 
averaged  and  displayed  as  a  frequency  by  a  Monsanto  Programmable  Counter-Timer 
Model  HOB. 

While  the  spin-up  rig  was  being  brought  up  to  the  desired  speed,  an 
air  jet(subsequently  referred  to  as  the  holding  jet)  from  a  nozzle  mounted 
rigidly  on  the  pedestal  at  the  end  opposite  the  bucket  wheel  held  the  obturator 
fixed  to  the  spinninq  rig.  The  holding  jet  passed  through  a  hole  in  the 
plunger  of  a  Rocker  Solenoid,  R.S.  No.  10-207.  This  hole  was  aligned  with 
the  hole  through  the  "projectile"  to  allow  the  jet  to  impinge  upon  the  obturator 
(Fig.  17).  When  the  appropriate  switch  was  activated,  standard  115VAC  was  applied 
to  the  Rocker  solenoid  and  removed  from  a  normally-closed  ASCO  Solenoid  Valve 
in  the  air  supply  line  to  the  holding- jet.  The  plunger  retracted  from  the 
position  shown  in  Fig.  17  and  the  solenoid  valve  closed.  The  retraction  of  the 
plunger  performed  three  functions.  It  first  caused  a  pulse  to  be  generated 
by  a  proximitor  mounted  next  to  the  nozzle.  This  pulse  started  the  timer 
function  of  a  second  Model  110B  Counter-Timer.  The  plunger,  by  misaligning  its 
hole  with  that  of  the  nozzle,  blocked  the  holding  jet.  This  removed  the  re¬ 
straining  influence  of  the  holding  jet  on  the  obturator  and  prevented  any  air 


remaining  in  the  supply  line  from  impinging  on  the  obturator.  And  lastly, 
in  the  fully  retracted  position,  the  machined  and  polished  end  of  the  plung¬ 
er  was  positioned  opposite  the  axis  of  the  apparatus  (Fig.  18).  When  the 
obturator,  now  free  to  move  relative  to  the  "projectile",  nutated  through 
enough  of  an  angle,  the  beam  of  a  SPECTRA  PHYSICS  MODEL  132  LASER,  MODEL  NO. 
3187,  passed  through  the  projectile  and  was  reflected  by  the  polished  end  of 
the  plunger  onto  a  light  sensitive  diode  (Fig.  19).  The  pulse  emitted  by  the 
diode  triggered  the  stop  channel  on  the  timer  and  the  elapsed  time  was  dis¬ 
played.  Figure  20  shows  the  entire  experimental  set-up  with  air  flasks  in  the 
background. 

EXPERIMENTAL  PROCEDURE 

The  obturator  was  positioned  in  the  spin-up  rig  in  such  a  way  that  the 

laser  beam  was  reflected  off  the  polished  spot  back  to  a  target  mounted  on  the 

laser  (Figs.  20  and  21).  The  center  of  the  target  is  the  location  of  the  beam 

and  the  obturator  was  adjusted  to  place  the  reflected  spot  on  the  center  of 

the  target.  Thus  the  orientation  of  the  z-axis  through  the  obturator  was 

known  relative  to  the  z-axis  (laser  beam)  through  the  projectile.  This  is  the 

initial  value  9. 

o 

Once  the  alignment  had  been  checked  and  the  obturator  positioned  and  held 
at  the  known  9q,  flow  to  the  turbine  nozzles  was  controlled  in  such  a  way  that 
the  desired  spin  rate  was  achieved.  The  turbine  flow  regulator  was  then  adjust 
ed  so  as  to  maintain  the  spin  rate  at  ±5  Hz  on  the  digital  display  or  approxi¬ 
mately  ±1.75  rad/sec.  Once  the  desired  rate  was  achieved  and  noted,  and  the 
timer  checked  and  reset  if  necessary,  the  solenoid  switch  (Fig.  17)  was  acti¬ 
vated.  The  obturator  was  released  and  began  to  move  relative  to  the  spin-up 
rig.  When  the  angle  9  reached  the  value  at  which  the  laser  beam  could  pass 
through  the  obturator  hole  [9=9*=  sin-1(r/R)],  the  beam  struck  the  polished 
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reflected  laser  beam  from  the  polished  end  of 
solenoid  plunger. 


plunger  and  was  reflected  into  the  light  sensitive  diode  housing.  The 
reflected  beam  striking  the  diode  caused  the  diode  to  emit  a  pulse  which 
turned  off  the  timer  gate.  The  elapsed  time,  from  plunger  retraction  to 
release  to  the  obturator  nutating  to  9*,  was  displayed  on  the  timer  dis¬ 
play  to  the  0.0001  second.  This  time,  t*,  was  recorded. 

EXPERIMENTAL  DETERMINATION  OF  THE  COEFFICIENT  OF  SLIDING  FRICTION,  y$. 

From  Eq.  (25)  it  will  be  noted  that  for  large  values  of  s&)pot  the  time 
elapsed  for  nutation  from  0Q  to  0*  is  given  by: 

t*  =  -[(eo  -  9*)  +  D]/(Rupo)  +  2?/a>po 

Without  the  influence  of  fluid  shear  (which  is  negligible  in  this  case), 

Eqs.  (34),  (35)  and  (39)  lead  to: 

t*  -*■  irAopoE(90  -  9*)  -  (A/2)sin28o]/(2ysWRsin0o) 

as  2sAopo-  0 

At  relatively  high  spin  rates,  therefore,  the  elapsed  time  t*  is  expected 
to  vary  directly  with  ajpQ.  Such  being  the  case,  the  coefficient  of  sliding 
friction,  may  be  estimated  from  the  slope  of  data  for  t*  vs  9*  at  high 
values  of  o> 

po 

DISCUSSION  OF  THE  EXPERIMENTAL  RESULTS 

Figure  22  shows  the  results  of  experiments  in  which  the  angular  velocity 
was  varied  within  the  range  of  capability  of  the  experimental  apparatus. 

Two  balls  were  tested,  each  with  different  hole  sizes  (see  Table  I)  but 
otherwise  essentially  identical.  Also  shown  in  Fig.  22  are  the  theoretical 
predictions  given  by  the  exact  solution  and  by  the  asymptotic  form  of  the 
linear  approximation,  Eq.  (41).  This  latter  expression  has  been  matched  to 
the  data  at  high  spin  rates  in  order  to  obtain  empirical  estimates  of  the 
coefficient  of  sliding  friction,  us- 
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Figure  22.  Presentation  of  experimental  results  showing  comparisons 
with  exact  and  approximate  solutions.  Gravity  load  only. 


The  experiments  verify  the  theoretical  prediction  of  the  high-speed 
asymptote  and  the  existence  of  a  minimum  in  the  elapsed  time  for  ball -open¬ 
ing,  with  respect  to  projectile  spin  rate.  The  discrepancies  between  theory 
and  experiment  at  lower  spin  rates  (where  the  linear  approximation  is  not  ex¬ 
pected  to  be  valid)  are  thought  to  be  mainly  due  to  inadequacies  in  the  ex¬ 
perimental  apparatus.  At  these  low  spin  rates  there  is  a  prolonged  period 
of  “hover"  of  the  ball  near  its  initial  position,  followed  by  a  rapid  de¬ 
crease  in  tilt  angle  (see  Fig.  8).  During  these  initial  instants  of  time 
any  external  disturbances  will  tend  to  cause  a  premature  conclusion  of  the 
hovering  period  and  it  is  this  effect  that  is  thought  to  have  led  to  the  low 
measured  values  of  t*  at  the  lower  values  of  ajp0.  The  effect  is  more  pro¬ 
nounced  in  the  Mod  I  ball  which,  due  to  its  small  value  of  X,  is  considerably 
less  prone  to  align  itself,  particularly  at  the  lower  values  of  cavity  spin 
rate. 

The  experiments  were  alsc  limited  in  precision  due  to  the  interaction 
of  the  holding  jet  with  the  solenoid-actuated  plunger.  This  difficulty 
contributed  to  the  scatter  of  data  and  limited  the  maximum  usuable  value  of 
cupo  because  of  the  high  jet  pressures  required  to  hold  the  ball  in  position 
at  high  spin  rates. 

The  empirically  determined  value  of  the  coefficient  of  sliding  friction, 
0.30  for  the  Mod  I  ball,  appears  to  be  reasonable  in  comparison  with  the  value 
of  0.35  reported  for  laminated  plastic  on  steel  [15].  The  larger  value  of  us  = 
0.51  for  the  Mod  0  ball  is  thought  to  be  due  to  the  influence  of  the  relatively 
large  hole  which,  in  the  high-speed  rotating  motion,  will  have  an  effect 
similar  to  that  of  roughness.  Within  the  precision  of  the  experiments  there 
was  no  detectable  influence  of  speed  upon  us  and  the  assumption  of  a  constant 
value  for  this  parameter  appears  to  have  been  justified,  at  least  for  the 
higher  values  of  <up0. 
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SENSITIVITY  TESTS  AND  DESIGN  GUIDES  FOR 
AERODYNAMIC  LOADING 


Figure  11  shows  a  comparison  of  the  linear  approximation  with  the 
exact  solution  for  the  motion  of  the  baseline  model  (Table  2)  at  the 
nominal  spin  rate  of  115,000  RPM.  Two  initial  ball  tilt  angles  are  also 
shown  and  it  may  be  seen  that  the  linear  approximation  yields  results  that 
are  quite  acceptable  for  initial  angles  near  90°.  The  values  of  M,  R,  and 
D  are  given  on  the  figure  for  reference  purposes.  The  usefulness  and  limita¬ 
tions  of  the  approximate  method  are  also  illustrated  in  Fig.  23  where  the 
initial  opening  time  is  given  as  a  function  of  projectile  spin  rate.  As 
in  the  experimental  program  this  time,  designated  as  t*,  is  the  time  required 
for  a  straight-through  path  to  open  along  the  projectile  centerline.  The 
ball  tilt  angle  at  this  instant  is  given  by 
9*  =  sin_1(r/R) 

and  the  corresponding  time  is  a  useful  reference  quantity  for  comparison  of 
various  obturator  designs.  (In  addition,  the  aerodynamic  force  model  used 
in  this  analysis  becomes  highly  suspect  at  values  of  9  below  9*  since  the 
simple  shock  structure  assumed  here  will  be  invalid  if  there  is  any  signifi¬ 
cant  flow  through  the  projectile). 

Figure  23  indicates  that  the  linear  approximation  approaches  the  exact 
solution  at  large  values  of  a)pQ.  This  is  to  be  expected  since  M  decreases 
quadratically  with  ojpp.  The  asymptotic  behavior  of  t*  at  high  a>pQ  is  a 
useful  design  insight  and  is  given  by 
t:  =  *opo[(90  -  9*)  +  D]/fisine0 

and  for  the  baseline  case  this  reduces  to  l*  a  4.0  x  10'4aj  ms 

00  po 

(<upo  in  s"1).  The  sharp  minimum  in  t*  illustrated  in  Fig.  23  suggests  some 


caution  in  the  design  of  these  devices  and  also  provides  encouragement  re¬ 
garding  the  existence  of  criteria  leading  to  an  optimum  design.  The  low 
spin -rate  behavior  shown  in  Fig.  23  is  outside  the  linear  range  of  operation 
and  is  due  to  the  delaying  mechanism  of  sliding  friction  at  relatively  large 
values  of  M  (large  friction,  low  inertial  imbalance  due  to  a  small  hole  in 
the  ball,  and/or  low  projectile  spin  rate). 

Under  these  conditions  the  ball  "hovers"  at  its  initial  value  and  con¬ 
siderable  time  elapses  before  a  ball  orientation  is  achieved  in  which  friction 
is  in  a  direction  to  aid  the  alignment  of  the  hole  with  the  projectile  axis. 
This  behavior  has  been  observed  in  previous  solutions  (see  Fig.  8),  and  in 
the  experimental  program  (see  Fig.  22).  Design  guidelines  necessary  to  avoid 
the  hovering  mode  are  discussed  below,  following  a  discussion  of  the  results 
of  the  sensitivity  analysis. 

SENSITIVITY  ANALYSIS 

The  design  variables  affecting  the  obturator  performance  are  those  con¬ 
tained  in  the  normalized  torque  determined  by  Eqs.  (18)  -  (21).  The  influence 
of  the  time-varying  relative  angular  velocity  and  moment  arms  in  these  ex¬ 
pressions  may  be  qualitatively  evaluated  by  referring  to  the  linear  approxi¬ 
mation  so  that  the  relevant  design  variables  are  [Eq.  40]  X,  ,  F=,  8,  9  . 

R,  and  <i)pQ.  Since  X  and  8  are  geometric  parameters  depending  upon  the  radii 
of  the  ball  and  its  hole;  and  the  minor  moment  of  inertia,  A,  depends  upon 
these  as  well  as  the  density  of  the  ball  material,  this  list  may  also  be 
written:  r,  R,  p  y  ,  M  ,  P  ,  and  co  Here  the  aerodynamic  force  has  been 

p  50000  no 

replaced  by  means  of  Eq.  (21)  and  the  projectile  hole  is  assumed  to  be  of  the 
same  size  as  that  of  the  ball.  In  considering  the  coefficient  y$,  it  must  be 
remembered  that  it  depends  upon  the  two  materials  (ball  and  projectile)  and 
is  therefore  dependent  upon  p.  We  have  studied  independent  variation  in  y$, 
however,  since  its  value  is  subject  to  considerable  uncertainty.  Experiments 
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have  shown  that  the  value  of  ps  should  be  approximately  50%  larger  than 
that  published  for  dry  friction  between  smooth  surfaces.  Table  2  indicates 
the  baseline  value  for  these  design  variables  as  well  as  the  maxima  and  minima 
that  determines  the  range  of  sensitivity  analysis.  In  each  case  (except 
“po’  whlc^  ^as  136611  Previously  discussed)  the  values  of  the  variables  were 
chosen  such  that  equal  increments  were  provided  and  the  baseline  condition 
was  included  in  the  range.  Single  parameter  variations  are  shown  here  al¬ 
though  work  is  continuing  using  the  techniques  of  numerical  optimization. 

The  computer  data  are  presented  as  9  vs  t  and  the  interested  reader 
may  construct  cross-plots  such  as  that  of  Fig.  23  to  evaluate  the  effects 
of  these  design  variables  upon  the  time  required  to  reach  a  given  angle,  such 
as  9*.  Figures  24-29  illustrate  the  relative  importance  of  the  selected 
design  variables.  It  is  important  to  note  that  in  these  figures  the  range  of 
system  performance  has  been  restricted  to  the  quasi -linear  range  previously 
discussed  where  the  ball  response  is  typically  of  the  type  shown  in  Fig.  8 
and  above  the  minimum  in  Fig.  23.  Outside  of  this  range  (in  the  "hovering" 
state)  the  trends  illustrated  may  be  reversed.  Attempts  to  reduce  ball  re¬ 
sponse  time  by  the  design  improvements  suggested  by  Figs.  24-29  may  push  the 
performance  beyond  the  minimum  illustrated  in  Fig.  23,  with  resulting  "negative 
improvements".  In  many  cases  the  implications  of  Figs.  24-29  are  self-evident 
but  to  assist  in  physical  interpretation  a  few  explanatory  remarks  are  offered 
below. 

Geometry  effects. 

Both  the  ball  radius,  R,  and  the  radius  of  the  ball  hole,  r,  were  varied 
independently.  Figure  24  shows  the  effects  of  increasing  the  ball  radius  while 
maintaining  a  constant  value  of  r/R.  The  variation  in  this  case  is  essentially 
one  of  scale-up  which  covers  a  range  of  projectile  diameters  from  20mm  to  100mm. 
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It  may  be  seen  that  larger  systems  respond  more  slowly.  Increases  in  applied 
torques  are  not  sufficient  to  offset  the  effects  of  the  fifth-power  relation¬ 
ship  between  ball  inertia  and  radius.  Figure  24  indicates  that  opening  of  the 
ball  obturator  in  large  projectiles  may  be  unacceptably  slow  unless  compensat¬ 
ing  design  actions  are  taken.  As  will  be  seen  in  subsequent  figures,  such 
actions  might  include  increases  in  ball  hole  radius  and  reduction  in  the  den¬ 
sity  of  the  ball  material. 

Increasing  the  size  of  the  hole  in  the  ball,  with  all  other  parameters 
held  constant,  reduces  the  inertia  of  the  ball  and  increases  the  imbalance 
between  the  moments  of  inertia  about  the  two  principal  axes.  Both  effects 
are  conducive  to  rapid  ball  response  and  this  is  illustrated  in  Fig.  25. 

In  addition,  the  initial  opening  angle,  9*,  is  directly  increased  by  an 
increase  in  r/R  with  a  corresponding  decrease  in  the  elapsed  time. 

Material  effects. 

Though  the  ball-obturator  must  be  of  sufficient  strength  to  pluq  the 
projectile  during  launch,  it  does  not  contribute  to  the  structural  integri¬ 
ty  of  the  projectile  in  flight.  The  ball  material  may  therefore  be  open 
to  some  design  variation  and  this  possibility  is  investigated  in  Fig. 26 
where  it  is  seen  that  the  time  required  to  reach  a  given  ball  angle  is 
approximately  proportional  to  the  density  of  the  ball  material,  p  A 
change  in  ball  material  from  Carbon  steel  (p  ^7, 750  kg/m3)  to  an  Aluminum 
alloy  (p  ^3,100  kg/m3)  would  lead  to  a  halving  of  the  response  time.  This 
effect  also  leads  to  an  opportunity  of  linearly  tailoring  the  ball  response 
time  by  means  of  its  density--a  feature  that  may  be  particularly  useful  in 
applications  of  the  system  as  a  hydromechanical  switch. 

Coefficient  of  sliding  friction. 

As  has  been  mentioned,  the  correct  value  of  is  subject  to  consider¬ 
able  uncertainty.  Figure  7  shows  the  impact  of  an  uncertainty  range  of 

about  ±20%  in  y$  which,  in  terms  of  t*,  is  about  ±15%.  In  applications  for 
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Figure  27.  Effect  of  coefficient  of  sliding  friction  on  ball  response. 
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which  margins  of  error  such  as  this  are  unacceptable,  tests  of  prototype 
configurations  may  be  necessary  to  determine  ys>  The  range  tested  in 
Fig.  27  represents  a  best  estimate  for  dry  friction  of  steel  on  steel 
for  the  range  of  r/R  covered  in  Fig.  25.  Systems  using  dissimilar  materi¬ 
als  will  require  appropriate  adjustments  in  yg.  In  addition,  it  may  be 
seen  that  some  roughening  of  the  Interfacial  surfaces  may  not  be  detri¬ 
mental  to  system  performance  while,  on  the  other  hand,  intential  polish¬ 
ing  of  these  surfaces  may  lead  to  slower  ball  action. 

Launch  conditions. 

Projectile  launch  Mach  number  and  altitude  have  been  separately  in¬ 
vestigated  with  the  results  shown  In  Figs.  28  and  29.  The  increase  in 
wave  drag  associated  with  higher  Mach  numbers  are  actually  beneficial  in 
terms  of  ball  response,  because  of  the  corresponding  Increases  in  the 
external  moments  that  aid  the  motion  of  the  ball  within  the  projectile. 

On  the  other  hand,  increases  in  launch  altitude  lead  to  Increases  in 
ball  response  times.  In  proceeding  from  sea-level  to  a  flight-level 
pressure  of  0.325  atm  (about  30,000  ft  in  the  standard  atmosphere),  for 
example,  a  threefold  increase  In  t*  is  predicted  for  the  baseline  design. 

DESIGN  GUIDES 

The  linear  approximation  previously  developed  represents  a  signifi¬ 
cant  simplification  to  the  computational  tasks  required  for  the  exact 
solutions.  As  such,  this  simplified  approach  is  recommended  for  design 
and  "first-cut"  calculations.  As  has  been  seen,  however,  there  is  a 
minimum  In  elapsed  time  (see  Fig.  23,  for  example)  which  represents  a 
demarkation  of  the  limit  of  validity  of  the  linear  approximation  as  well 
as  a  point  of  optimum  design  for  the  ball  obturator.  A  means  for  esti¬ 
mating  this  optimal  point  is  developed  below. 
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Under  the  initial  conditions  of  the  problem  at  hand  (tj/=l,  1=0,  0=0, 
e=eQ)  the  equations  of  motion  [Eqs.  (8)  -  (10)]  may  be  written: 

$  *  e0cote0 

$  *  -6  cos26  /sine 
0  0  0 

9  .  _M6o  -  (X/2)sin20o 

Here  the  term  M6q  is  retained  since  It  may  be  significant  for  large  values 

of  M  even  though  9q  is  initially  zero.  In  addition,  in  the  initial  instants 

of  motion,  the  relative  motion  of  the  two  bodies  is  largely  due  to  nutation. 

That  is,  oj  =  -8,  and 
r  o 

M0  =  -M/Aup02 

Thus  an  approximation  to  the  Initial  acceleration  0  is  given  by 
0  a  ( M/ Au)p0 2 )  +  0  =  constant 

In  cases  in  which  hovering  occurs  this  expression  will  be  approximately 
zero  and  an  Indication  of  a  hovering  design  is  therefore  given  by 

po 

and  the  design  guide  to  avoid  hovering  Is  given  by 

_ R  <  ! 

In  the  initial  Instants,  the  applied  moment  due  to  aerodynamic  loading 
may  be  conservatively  estimated  as  R$a  ■  vsFaRcos8  so  that  the  design 
criterion  Is  given  as 

2  u  F  RcosS 

H  =  - -  <  1  (42) 

XAw  2s1n20_ 
po  o 

where,  as  usual,  8  *  s1n-1(r/R).  Note  that  In  this  expression  as  X  -*•  0 
or  0  -►  ir/2  the  factor  H  grows  beyond  bound.  The  first  situation  occurs 


If  there  is  no  hole  in  the  ball  (A  *  C)  and  the  second  case  represents  the 
stable  initial  condition.  Equation  (42)  indicates  a  hovering  condition  in 
both  cases,  in  accordance  with  expectations. 

As  an  example,  consider  the  baseline  case  (Table  II).  For  this 
design  and  Mach  number,  F.  ~  93  N  (about  21  lb).  Inserting  the  other 

Q 

constants  from  Table  II  we  have 
H  -  0.735 

and  the  prediction  is  that  hovering  is  avoided.  If  the  same  design  is  con¬ 
sidered  but  Eq.  (42)  is  solved  for  u>^Q  in  the  critical  case  (H  =  1)  we 
have  a)pQ(critical )  =  10,000  s'1  or  about  98,000  rpm.  Comparison  of  these 
calculations  with  the  results  indicated  in  Figs.  8  and  23  shows  a  very 
good  agreement.  In  particular,  Fig.  23  shows  that  if  the  equality  is 
used  in  Eq.  (42)  a  very  good  estimate  of  the  optimum  design  is  obtained. 
Thus,  for  a  given  worst-case  launch  condition  (maximum  F_  due  to  maximum 

a 

Mach  number  and  minimum  altitude)  and  initial  projectile  spin  rate  (low¬ 
est  value  for  worst  case)  the  optimum  ball  design  (X,  A,  8,  y  ,  0Q)  may 
be  deduced  from  Eq.  (42).  It  would  be  wise,  of  course,  to  apply  some 
safety  factor  to  this  optimum  design  and,  in  any  case,  a  violation  of 
the  inequality  in  Eq.  (42)  would  lead  to  hovering  -  a  failure  in  the 
operation  of  the  ball -obturator.  To  repeat,  if  the  inequality  is  satis¬ 
fied  in  Eq.  (42)  the  linear  approximation  [Eqs.  (25)  and  (40)]  may  be 
safely  applied. 
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CONCLUSIONS 


The  gyrodynamic  motion  of  a  ball  obturator  within  a  spinning  tubular 
projectile  has  been  investigated  from  both  theoretical  and  experimental 
points  of  view.  The  theoretical  developments  have  included  both  exact  and 
approximate  formulations  and  comparison  of  these  two  theories  has  led  to 
the  definition  of  the  range  of  validity  of  the  approximate  method. 

The  experimental  program  was  conducted  under  conditions  in  which  the 
main  source  of  external  moment  was  sliding  friction  due  to  the  weight  of 
the  ball.  The  results  of  these  tests  have  substantiated  the  validity  of 
both  the  exact  and  approximate  theories  and  have  led  to  several  insights  re¬ 
garding  the  motion.  In  addition,  a  semi -empirical  method  of  measuring  the 
coefficient  of  sliding  friction  evolved  from  the  experimental  program. 

Two  distinct  regions  of  motion  have  been  identified.  The  first  of 
these  -  the  hovering  mode  -  is  characterized  by  a  metastable  initial  ori¬ 
entation  for  prolonged  periods  of  time.  This  undesirable  behavior  may 
be  avoided  if  design  guides  developed  within  this  study  are  followed. 

The  second  region  of  motion  occurs  when  the  hovering  region  is  avoid¬ 
ed  and  in  this  mode,  the  response  of  the  ball  is  approximately  linear.  The 
analysis  has  shown  that  minimum  ball -opening  times  may  be  expected  for  ball 
designs  and  projectile  launch  conditions  that  are  near  (but  not  beyond)  the 
onset  of  hovering  motion. 

The  analytical  model  has  been  used  to  conduct  sensitivity  tests  to  de¬ 
termine  the  relative  influence  of  design  parameters  for  ball  designs  that  are 
in  the  proper  (non-hovering)  range.  In  general,  ball  response  is  more  rapid 
for  larger  hole  sizes,  lighter  material,  greater  sliding  friction,  higher 
launch  Mach  number,  and  lower  launch  altitudes.  Small  obturating  balls  will 
respond  more  quickly  than  large  ones,  all  else  being  equal.  The  theoretical 


predictions  depend  directly  upon  the  nature  and  magnitude  of  the  aero¬ 
dynamic  load  acting  upon  the  ball  in  flight.  Uncertainties  regarding  the 
theoretical  description  of  this  load  are  the  main  source  of  doubt  re¬ 
garding  the  precision  of  the  predictive  models  presented  herein. 
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LIST  OF  SYMBOLS 


A  Mass  moment  of  inertia  about  minor  axis 

C  Mass  moment  of  inertia  about  major  axis  (axis  of  hole  in  ball) 

D  Step  input  associated  with  linear  model 
£  Unit  vector 

F  Force 

h  Gap  width  separating  ball  and  cavity  (fluid  friction  model) 

H  Design  criterion,  see  Eq.  (42) 

k  Ratio  of  specific  heats 

M  External  moment  (also  used  to  denote  Mach  number) 

P  Pressure 

r  Radius  of  hole  in  ball 

R  Radius  of  ball 

J?  Slope  of  ramp  input  associated  with  linear  model 

r  Moment  arm 

m 

t  Time 

Greek  Letters 

0  Angle  used  in  locating  contact  points  (  =  9*  =  sin_1r/R  ) 

Y  Angle  used  in  locating  contact  points  (see  Fig.  5) 

e  Dimensionless  gap  width  (  =  h/R  ) 

S  Damping  coefficient  associated  with  linear  model 

n  Angle  used  in  locating  contact  points  (see  Fig.  5) 

X  (C  -  A)/A 

y  Molecular  viscosity  of  fluid  in  gap 
p  Density  of  ball  material 
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LIST  OF  SYMBOLS 


$> 

<M»e 

U) 


1,2,3 

a 

f 

g 

n 

o 

r 

s 

x,y,z 


Greek  Letters  (continued) 

Phase  angle  associated  with  linear  approximation  (see  Eq.  25) 

Euler  angles 
Angular  velocity 
Spin  rate  of  projectile 

Subscripts 

Refer  to  body  fixed  axes 
Aerodynamic  quantity 
Fluid  quantity 
Refers  to  gravity  load 
Normal  component 

Initial  quantity  (also  refers  to  solutions  of  torque-free  motion) 

Denotes  relative  angular  velocity 
Associated  with  sliding  friction 
Refer  to  inertial  (projectile-fixed)  axes 
Vector  quantity 
Conditions  in  free  stream 

Superscripts 

Dimensional  quantity  -  absence  infers  non-dimensional i zed  time-dependent 
quantity.  E.g.  0  =  S((upo)  • 

Refers  to  condition  of  first  opening  of  line-of-sight  along  projectile 
axis  and  through  hole  in  ball  obturator 
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